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Exact solutions of the Dirac equation in external electromagnetic background fields are very help¬ 
ful for understanding non-perturbative phenomena in quantnm electrodynamics (QED). However, 
for the limited set of known solutions, the field often depends on one coordinate only, which could be 
the time t, a spatial coordinate such as x or r, or a light-cone coordinate such as ct — x. By swapping 
the roles of known and unknown quantities in the Dirac equation, we are able to generate families of 
solutions of the Dirac equation in the presence of genuinely space-time dependent electromagnetic 
fields in 1 -I- 1 and 2-1-1 dimensions. 
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I. INTRODUCTION 

Quantum electrodynamics (QED) as the theory of 
charged particles interacting with electromagnetic helds 
is well understood in the context of standard perturba¬ 
tion theory and can describe several intriguing phenom¬ 
ena of nature. However, QED contains other fascinat¬ 
ing effects that cannot be explained using perturbative 
methods. Such non-perturbative effects can arise when 
the electromagnetic field is so strong that it cannot be 
treated as a perturbation. In order to understand these 
phenomena, it is often useful to study the behaviour of 
exact solutions in those external background fields. 

Unfortunately, although the Dirac equation was first 
formulated more than eighty years ago the set of 
known exact solutions is still quite limited (see, e.g., [2] 
for a review). Apart from the Coulomb field oc 1/r^ 
Ellll, exact solutions are known for a constant electric 
field and a Sauter profile in space oc 1/cosh^(fcx) or time 
oc 1/ cosh^ (wt), for example. The latter are relevant for 
the non-perturbative Sauter-Schwinger effect mm corre¬ 
sponding to electron-positron pair creation from vacuum 
via tunnelling. In contrast to electron-positron pair cre¬ 
ation in the perturbative (multi-photon) regime which 
has been observed at SLAG [5], this non-perturbative 
prediction of quantum field theory has not been conclu¬ 
sively experimentally verified yet. However, there are 
several experimental initiatives which might be able to 
eventually reach the ultra-strong field regime necessary 
for observing this striking effect [5]. 

Furthermore, exact solutions are known for a constant 
magnetic field (relativistic Landau levels, see [TUI ITT] ) 
and plane waves, where the fields depend on one of the 
light-cone coordinates such as ct — x (Volkov solutions, 
see e.g. [TUHT5] ). These (transverse) helds do not induce 
pair creation from vacuum. 

Nevertheless, in all these cases, the helds depend on 
one coordinate only (such as r, x, t, or ct — x). As a 
result of this high degree of symmetry, the set of partial 


differential equations can be reduced to an ordinary dif¬ 
ferential equation, which greatly simplihes the analysis. 
An analogous limitation applies to our theoretical un¬ 
derstanding of the Sauter-Schwinger ehect. Even though 
there are many results for helds which depend on one 
coordinate only, we are just beginning to understand the 
impact of the interplay between spatial and temporal de¬ 
pendencies, see, e.g., [MU]. 

In the following, we develop a method which allows us 
to obtain solutions of the Dirac equation for genuinely 
space-time dependent helds. To this end, we pursue a 
different approach by assuming that we already know a 
solution to the Dirac equation. We then calculate the 
vector potential corresponding to the given solution 
from the Dirac equation. This is feasible as the Dirac 
equation does not contain any derivatives of the vector 
potential. More generally speaking, we write down a so¬ 
lution to a partial diherential equation and then try to 
hnd a physical problem associated with the solution - a 
concept also well known in the held of huid dynamics, 
see, for example, 


II. LIGHT CONE COORDINATES 


Let us start with the most simple and yet non-trivial 
case - the Dirac equation in 1-1-1 dimensions. For the 
following derivation, it is convenient to transform to light 
cone coordinates x± dehned as (?i = c = 1) 


t + X 


t — X 


( 1 ) 


Since perhaps not all readers will be familiar with the 
form of the subsequent expressions in light cone coordi¬ 
nates, let us insert a brief reminder. The Jacobian matrix 
of the coordinate transformation between Cartesian and 
light cone coordinates 




d{x+,x_) 
d (t, x) 


1 
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( 2 ) 
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yields the transformation laws for tensors such as the 
partial derivatives 


d' = {J 


-l\iy 


Y d 


1 

71 


(dt + dx 

\dt-dx 




( 3 ) 

In 1+1 dimensions, the electromagnetic field strength 
tensor contains only one independent component, the 
electric field E{t,x) 


which thus reads in light cone coordinates 

• ( 5 ) 


Transforming the Cartesian Minkowski metric tensor 
ilfj.i' = diag(+l, —1) to light cone coordinates as well gives 

9pv = J) • (6) 

A possible choice of light cone gamma matrices satisfying 
the Clifford algebra’s anti-commutation relation 

{7^7"} = 25'^", (7) 

therefore is 

= (o f= (^^ o) • 

Note that in 1+1 and 2+1 dimensions, the Clifford alge¬ 
bra can be satished with 2 x 2-matrices. 


III. INVERSE APPROACH 


For arbitrary these expressions are not necessarily real. 
Therefore, we require the imaginary parts of qA^ and 
qA_ to vanish, giving two conditions which we use to 
eliminate two real degrees of freedom of the spinor ip. 
Using the polar representation for the spinor components 
'ipk = these conditions can be written as 


Tfl 

r 2 -^ rir 2 sin [ipi - tpa) = 0, 

rn (12) 

Ti d_ri + —= rir 2 sm {(pi - (^ 2 ) = 0. 
v2 

Adding the two equations gives 


d_rl = -d+rl 


m 


d_ri = r 2 sin [p 2 - + 1 ) • 


(13a) 

(13b) 


The first equation (13a) can be solved for r 2 by integrat¬ 
ing with respect to x+ 


r2 = 



d_ rl dx+ , 


(14) 


where c{x-) is an integration constant that may still de¬ 
pend on X-. The remaining equation (13b) determines 
the phase difference +2 ~ +1 


+2 — +1 = arcsm 


\/2d_ri 
m r 2 


(15) 


where we could also use other branches of the arcsin- 
function such as Ap = +2 ~ +1 —>■ tt — Ap, leading to dif¬ 
ferent solutions in general - see the remark after Eq. (17). 
Using the abbreviation 


The Dirac equation, minimally coupled to the electro¬ 
magnetic potential via the charge q 

(* 7 '" [5^ +iqA^,] -m)ip = 0, (9) 

assumes the following simple form in terms of the light 
cone gamma matrices (1^ 


—m 


i\/2 [i9_ +iqA-] 


i\/2 [a_,_ +igA+]\ (ipi 


—m 


) 


= 0. (10) 


Traditionally, the Dirac equation is treated as a partial 
differential equation. A solution ip for a specific potential 
A^ is typically calculated by reducing the Dirac equation 
to an ordinary differential equation. In our approach, we 
assume that we know a specific spinor ip = {ipi,ip 2 )'^ 
which is a solution to the Dirac equation and calculate 
the corresponding potential. Thus, we solve ( [l0| for the 
components of 


qA+ = i + 
qA- = i 


d+ ip2 


'4’2 

d ipi 




m ipi 

V21P2'' 

m ip 2 
y/2^Pi' 


( 11 ) 


= l/c- / d_r^dx+ - ^{d_ry 


(16) 


and Eqs. (14) and (15), we can calculate the form of the 
spinor ip 


Ip = 


±s ■ 


r 


d_ 


(17) 


where we have set r = ri and tp = tpi. Note that we 
find two different solutions with ±s corresponding to the 
different branches of the arcsin or square-root functions 


in Eqs. (15) and (16), respectively. 


Local gauge invariance allows us to eliminate the phase 
6*'^’ by applying a gauge transformation ip ip' = e~''‘'^ip, 
which adds a term p to qA^. The components of A^ 
using the spinor given in 0 finally are 


, m r \/2 i9 1 d_ r 

9^+=T^- +-^-, 

Y2 sms 

A ^ s 


( 18 ) 
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These are obviously real as long as r and s are real, too. 
The electric field corresponding to this potential accord¬ 
ing to © is 


E = d_A+-d^A_. (19) 

In summary, by choosing a real generating function 
r{x+,x-) and a real supplementary boundary value func¬ 
tion c(x-), we can generate arbitrary space-time depen¬ 
dent solutions 'i/’(x-i., X-) of the Dirac equation in the 
presence of an electromagnetic background which 

can also depend on space and time. 

Obviously, the associated electromagnetic field 
strength tensor in Eq. ([^ automatically satisfies the 
homogeneous Maxwell equations as it has been derived 
from a vector potential A^. If we demand that it also 
obeys the inhomogeneous Maxwell equations (/tq = 1) 


Thus, a gauge transformation tp ^ ip' = e with 


Pm 



m fr/s\ 

y/2 V«/v 


(24) 


can be used to set the potential components to zero and 
reveals that these solutions are plane wave solutions to 
the free Dirac equation of either positive or negative en¬ 
ergy. Transforming the light-cone momenta p± back to 
the usual Cartesian representation po = {p+ + P-) /V^, 
we find that the energy is given by po = ra{r/s -\- slr)l2. 
Thus a solution where both r and s are positive (or both 
negative) corresponds to a positive energy whereas dif¬ 
ferent signs of r and s yield a negative energy. 


B. Single pulses 


( 20 ) 


we have to specify the sources j'' accordingly. In 1 -|- 1 
dimensions we find 


P= -dxE 
J= d,E 


= ^{d_E-d^E), 
= ^{d_E + d^E), 


( 21 ) 


where p is the charge density and j is the current density. 
For non-trivial field profiles E{t^x), they will be non¬ 
zero in general. However, this is no surprise because the 
only vacuum solution of the Maxwell equations in 1 -|- 1 
dimensions is a constant electric field E = const. 


IV. SOLUTIONS 

In order to illustrate the approach presented in the 
previous section, let us discuss some exemplary solutions 
that can be found using this method, starting with the 
most simple ones. The expressions for the spinor and the 
potential components are significantly simplified if r is 
independent of a;_. 


In this subsection, we find solutions for arbitrary light 
cone fields E(x+) and E{x-), i.e., pulses moving along 
the light lines. Such solutions were found before using 
traditional methods as well [TiirB] . 


1. x+-dependent pulse 

Let us assume that the function r depends on x+ only 
while s = ±i/c = const 


iP = 


’(a^-s) 


(25) 


In this case, neither the spinor nor the vector potential 
depends on x_ which simplifies the expression for the 
electric field 


777 1 

qE = qd_ A_|_ —q d_^_ A_ = 


V2 r{x+)' 


(26) 


=0 


This is a first-order ordinary differential equation for 
r(x+) which can be integrated easily 


-(x+) = ' 


l + yl^Eq r E{i+)dx+ 
m s 


, (27) 


A. Plane waves 


Choosing r and s = ±^/c to be constant. 


iP 



= const, 


( 22 ) 


leads to a constant electromagnetic vector potential 


qA+ = 
qA_ = 


m r 

-- = const, 

72 5 
m s 

-- = const. 

V2r 


(23) 


with Tin = r(x+ —>■ —oo). Comparison with Sec. IV A re¬ 
veals that the pre-factor •\/2?'in/ (ms) in front of the above 
x+-integral over qE is just the inverse initial momen¬ 
tum 1/p™. As already discussed in [T3], the term in the 
square bracket in Eq. (271 vanishes and thus r diverges 
when this x+-integral over qE becomes large enough to 
compensate p™. Note that the light cone dispersion rela¬ 
tion p+p- = m^/2 shows that must diverge when p_ 
vanishes and vice versa. 

Now, let us recall that the phenomenon of electron- 
positron pair creation (such as in the Sauter-Schwinger 
effect) can be described by the situation where an ini¬ 
tial solution with positive energy transforms into a final 
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solution which contains contributions with negative en¬ 
ergies (or vice versa). Assuming that r becomes constant 
initially and finally, we find that pair creation can only 
occur if r{x+) changes its sign somewhere, i.e., if r{x+) 
vanishes or diverges at some point. If r(a;+) crosses zero, 
the electric field (26) diverges - whereas a diverging r{x+) 
precisely corresponds to the case discussed above, see also 
M- Thus, we find that we cannot describe particle cre¬ 
ation in this case without introducing some singularity 
(see the Appendix). 


2. X--dependent pulse 


In a similar way, we can derive solutions for electric 
fields only depending on X- by setting r = const and 
letting s(a;_) = ±-\/c(a;_) depend on a;_ 




(28) 


Thus, the electric field can be calculated as follows 

777 1 

qE = qd_A+-qd,A_ = - -=rd_ 


s[x_y 


(29) 


=0 


which is again a first-order ordinary differential equation 
for s(a:_). The solution is given by 


s(x_) = 


r~ E{i.)dx_ 


, (30) 


with Sin = s{x- —> —oo). In complete analogy, the same 
arguments as for an x_|_-dependent pulse apply in this 
case. 


C. Two pulses 


As a non-trivial extension of these two cases, we can 
combine the previous two solutions into a single spinor 




(31) 


where the two components are given by 


"{x+) =■ 


5(X_) = . 


f ^ E+iS:+)dx+ 
m Sin 7-00 

r~ E_{i_)dx. 


-\ -1 


m ri 


in J —oo 


-1 


(32) 


We may calculate the electric field using (18) and (19) 


E{x+,X-) = '-^^E_{x-) + ^^^E+{x+). (33) 



FIG. 1. Plot of r{x+,X-) as given in (341 with rin = 1, 
5 = 0 . 2 , and 7 = 1.2/m. 


Initially, we have E{x+,X-) = E_{x-) if+(x_|_) which 

corresponds to two independent pulses approaching each 
other from different directions. When these two pulses 
meet, however, this is no longer true - which shows that 
the mapping from '0 (i-e., r and s) to A^ is not linear. For 
late times, these pulses propagate again independently, 
but with modified amplitudes in general. 


D. Emerging pulses 


Another solution where the corresponding electric field 
consists of two pulses can be generated by setting 


'■(a:+,x_) = Tin -I- 




1 -I- e-^^+ -I- e-T'^- 


(34) 


For non-vanishing ^ and 7 > 0, the chosen r(x+, x-) will 
be constant almost everywhere except in the vicinity of 
the forward light cone (see figure [^. 


In this case, the expression for s according to (16 1 is 


not as simple as before because r is not independent of 
X-. Nevertheless, s can be calculated analytically, al¬ 
though the resulting expressions for s and the electric 
field qE are quite lengthy. Thus, we will only give a plot 
of the resulting electric field which shows the two pulses 
emerging from the origin and moving along the forward 
light lines (see figure [^. 


V. EXTENSION TO 2+1 DIMENSIONS 

The approach presented here can be extended to 2 + 1 
dimensional space-times as well. We use the Cartesian 
coordinate y in addition to the light cone coordinates x+ 
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is a solution of the Dirac equation with the potential 
components 


qA+ = 


m ri 

y /2 S 


1 5 ri d, u 


m r| 
a /2 ris 


qAy = —m 


u 

s 


1 ^2 dy r 2 u d_ ri 

a /2 riS ris 

udyTi 1 9_,_r| 

riS y/2 ris 


(41) 


In contrast to 1+1 dimensions, the electromagnetic field 
strength tensor contains three independent components, 
for example the two electric fields E^^y in x and y di¬ 
rection plus the perpendicular magnetic field B^- These 
components of the electromagnetic field can be calculated 
as follows 


FIG. 2. Plot of the electric field qE corresponding to the 
solution generated by r{x+,x-) given in (341 with rin = 1, 
^ = 0 . 2 , and 7 = XAjm. 


and X-. Thus, the metric tensor becomes 


=d_A+- dj^ , 

Ey = {d_ Ay — dy A_ + dj^Ay — dy A_|_) , 

Bz = ^ {d_ Ay — dy A- — d_f_Ay + dy A_|_) . 


/O 1 0 \ 

5 '^" = 1 0 0 . (35) 

VO 0 - 1 / 

In order to complete our set of gamma matrices from (|^ , 
we choose the third 7 -matrix according to 

(36) 

Thus the Dirac equation in 2 + 1 dimensions is given by 

/-TO - [dy +iqAy] iy/2 [9+ EiqA+\ \ /^ ^ 
\i'/^[d_+iqAJ\ -m+\dy+iqAy\) yl^2) 

(37) 

In complete analogy to section m we solve the Dirac 
equation for qA^ and qA_ and reduce the spinor’s num¬ 
ber of degrees of freedom by requiring the imaginary 
parts of the electromagnetic potential’s components to 
vanish. After some calculation, we are able to write the 
spinor and the electromagnetic potential in terms of three 
real functions ri{x+,x-,y), r 2 {x+,x-,y) and c{x+,x-). 
Explicitly, a spinor of the form 



We see that these expressions simplify significantly if ri 
and r 2 are independent of y. In that case, the electromag¬ 
netic field does only depend on the light cone coordinates 
as before and similar solutions as in the I +1 dimensional 
case can be found, e.g. one and two wavefronts. In fact, 
the solutions given in section [TV] are solutions to the 2 + 1 
dimensional Dirac equation as well but can be extended 
to also include a transverse electric and magnetic field 
component. 

To verify that our method reproduces known solutions, 
we insert the lowest Landau level solution 


/ = A/ exp 





(43) 


into our formalism, i.e. we set 

ri{x+,x-) = TV exp ^“^ 9 ^ 
r 2 {x+,x-) =ri{x+,x-), c = 0 , 

(44) 

where A/ is a normalization constant. Calculating the 
potential components gives 


x+ — X- 

~w~ 


-[ 2 N 


qB 


with 


s = ± 



(39) 


qA+ = qA_ = qAy = -qB ^^^ + ky, (45) 

so that the electromagnetic field is 


and 


y 2 ri 


2 {x+,X-)+ J {d_rj+d_^_rl) dy 


(40) 


Ex = Ey = 0, Bz = B, 
which is the expected result. 


( 46 ) 
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VI. CONCLUSIONS & OUTLOOK 


Appendix A: Perturbed solution 


We have developed an inverse approach for generating 
families of exact solutions of the Dirac equation in the 
presence of space-time dependent electromagnetic fields 
in 1 - 1-1 and 2-1-1 dimensions. Somewhat similar to op¬ 
timal control theory, we start with a suitable ansatz for 
the spinor ip and then derive the appropriate background 
field Afj^ which supports this solution. In 1-f 1 dimensions, 
we may choose a real generating function r(x+,x-) and 
a suitable real supplementary boundary value function 
c(x_) such that the radicand in Eq. (14) stays positive. 
In 2-I-I dimensions, we may choose two real generating 
functions ri{x+,x-,y) and r 2 {x+,X-,y) as well as one 
real boundary value function c{x+,x-). 


The solutions generated in this way may depend on 
space and time in a complicated manner - a situation 
which is quite difficult to treat with traditional meth¬ 
ods. As one possible application, our method could be 
used to solve steering problems such as: given an initial 
wave-packet tpim which electromagnetic field induces 
an evolution to a prescribed final wave-packet V'out? As 
another application, these exact solutions could be used 
as touchstones for already existing exact or approximate 
non-perturbative derivation techniques (e.g. the world¬ 
line instanton method [23] 1 or as starting point for new 
approximative methods, such as WKB [21] or lineariza¬ 
tion around a given background solution (see the Ap¬ 
pendix). 


The structure of the Dirac equation suggests that this 
general strategy can also be applied to 3-1-1 dimensions, 
where both the potential and the Dirac bi-spinor 
Ip have four components. Thus, for a given ip, we get 
four equations for the four components A^, which can be 
solved (except in singular cases). However, the four con¬ 
straints 5(A^) = 0 assume a form which is far more com¬ 
plicated than in 1-1-1 and 2-1-1 dimensions. This renders 
the identification of real generating functions which cor¬ 
respond to the remaining degrees of freedom rather cum¬ 
bersome. The analysis could be simplified by restricting 
the space-time dependence to l-fl and 2 - 1-1 dimensions, 
which should be the subject of further investigations. 
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To find solutions for electric fields that create electron- 
positron pairs (see, e.g., [25H28] 1. we use the ansatz 

r = a + P sin{mj), (Al) 


where the Bogoliubov coefficients a and P as well as the 
eikonal function 7 are slowly varying functions of the 
light cone coordinates. (We consider 1-1-1 dimensions for 
simplicity.) The main idea here is that a is an exact 
solution and P is used to slowly turn on an oscillating 
perturbation. The value of P then is related to the pair 
creation rate. 

However, the calculation of s and qE is rather com¬ 
plicated for arbitrary functions a, P, and 7 because s 
depends nonlinearly on r. Hence, as the perturbation 
should be small, we calculate the electric field only up to 
linear order in P 


E = -t -t o{p^), 


(A2) 


where qE^°‘^ is the unperturbed force of order /3° and 
qE^^'> is the first-order perturbation of order P^. Apart 
from this linearization, we assume that the mass m rep¬ 
resents the largest energy scale in the problem and thus 
we employ a large-m expansion on top of the approxima¬ 
tion in Eq. (A2). Expanding qE^^'> into powers of m and 
keeping only the highest-order term gives 




rrP{dj^ 7 )^( 5 - lY 


Sa d+1 

m a m Sa ^ \ , 1 

— 5 _ 7 H--=— 0+7 ' 

y/2Sa y/2 a ^ 




with the abbreviation 


= -qA 


(c) 


-I- 0{m ), 
(A3) 


= \ c — d dxA 


(A4) 


(gA^^ and qA'pp^ are the leading-order contributions to 
the vector potential.) Since a, P, and 7 are s upp osed to 
be slowly varying, the leading contribution (A31 would 
be rapidly oscillating due to the pre-factor cosimj) un¬ 
less the phase function 7 has a stationary point (see be¬ 
low) . Of course, such a rapidly oscillating force with a fre¬ 
quency of order m could well create pairs, but this process 
would be typically in the perturbative (multi-photon) 
regime. Here, we are interested in non-perturbative phe¬ 
nomena such as the Sauter-Schwinger effect and thus we 
demand that these rapidly oscillating contributions are 
absent - at least to leading order. Thus, we require the 
term of order rnP in qE^^'> to vanish. This is the case if 
S = m"f solves the eikonal equation 


1(0 


= (a^5 + gAV“’) [d_S + qA[ 


(A5) 
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Therefore, this condition can be used to fix 7 for a given 


a. Then, the leading order of is of order 



where we have used the eikonal equation (A5) to sim- some point. As those are proportional to 


plify some expressions. If we require this rapidly oscil¬ 
lating term to vanish as well, we get a linear first-order 
partial differential equation for /3. However, this linear 
equation does not have any source term. Therefore, a 
solution where j3 vanishes initially will not generate any 
pairs unless the coefficients of 9_|_ /3 and d_ /3 vanish at 


m a , m Sa „ 


(A7) 


we only obtain pair creation at this level of description 
if (A7| vanishes somewhere. According to the eikonal 


equation (A3), this in turn implies that 9,7 or 9_ 7 


has to vanish somewhere, i.e., that the phase function 
becomes stationary. 
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